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The Small Scatter in BH-Host Correlations & The Case for 
Self-Regulated BH Growth 
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ABSTRACT 

Supermassive black holes (BHs) obey tight scaling relations between their mass and their 
host galaxy properties such as total stellar mass, velocity dispersion, and potential well depth. 
This has led to the development of self-regulated models for BH growth, in which feedback 
from the central BH halts its own growth upon reaching a critical threshold. However, models 
have also been proposed in which feedback plays no role: so long as a fixed fraction of the host 
gas supply is accreted, relations like those observed can be reproduced. Here, we argue that the 
scatter in the observed BH-host correlations presents a demanding constraint on any model 
for these correlations, and that it favors self-regulated models of BH growth. We show that 
the scatter in the stellar mass fraction within a radius R in observed ellipticals and spheroids 
increases strongly at small R. At fixed total stellar mass (or host velocity dispersion), on very 
small scales near the BH radius of influence, there is an order-of-magnitude scatter in the 
amount of gas that must have entered and formed stars. In short, the BH appears to "know 
more" about the global host galaxy potential on large scales than the stars and gas supply 
on small scales. This is predicted in self-regulated models; however, models where there is 
no feedback would generically predict order-of-magnitude scatter in the BH-host correlations. 
Likewise, models in which the BH feedback in the "bright" mode does not regulate the growth 
of the BH itself, but sets the stellar mass of the galaxy by inducing star formation or blowing 
out a mass in gas much larger than the galaxy stellar mass, are difficult to reconcile with the 
scatter on small scales. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — quasars: general 
— cosmology: theory 



1 INTRODUCTION 

The tight empirical correlations between the masses of super- 
massive black holes (BHs) and the velocity dispersion, masses, 
and binding energy or potential well depth of their hosts demon- 
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). Understanding the physical origin and consequences of 



these correlations is critical for informing models of the co- 
formation of BHs and bulges, as well as theories which relate the 
evolution and statistics of BH formation and quasar activity to the 
remnant spheroid population. Likewise, the interpretation of obser- 
vations tracing the buildup of spheroid populations and associations 



* E-mail:phopkins @ astro.berkeley.edu 



between spheroids in formation, mergers, and quasar hosts depends 
on understanding the evolution of the BH-host correlations. 

Although the characteristic spatial and mass scales of the BH 
and tightly correlated host properties are wildly different, their 
characteristic energy and momentum scales are the same. That is, 
a few percent of the radiated luminosity or momentum from the 
BH growth is comparable to the binding energy/momentum of 
the galactic gas supply. Motivated by this comparison, attempts 
to explain these correlations have led to the development of a 
class of self-regulating feedback models, which argue that the en- 
ergy or momentum released from an accreting supermassive black 
hole, even if only a small fraction couples to the surrounding 
ISM, is sufficient to halt further accretion onto the black hole 
and drive away gas, self-regulating growth by shutting off the 
qua sar and potentially quenching star formation in t he host, (see, 
e.g.. lCiotti & OstrikeJl997ll200lllsiik"& Reesll998l : lMurrav et all 
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2005|). In these models, star forma- 



tion and inflows proceed rapidly before the final stages of BH 
growth, so that the BH grows in a relatively fixed background po- 
tential, which sets the critical BH accretion rate at which the BH 
halts its own subsequent growth[] Implicit in these models, the idea 
that the BH growth is fueling in sudden, violent events (mergers) 
that lead to strong bulge and BH growth, links the formation of 
both (it is less clear whether or not such models could succeed if 
growth occurred more slowly, in low accretion rate states, where 
there might be little bulge formation and/or feedback). 

The gas mass blown out by the BH may be an order of magni- 
tude larger than the BH mass itself, but is still a small fraction of the 
total galaxy mass; it is the mass "left over" from a central starburst. 
This picture is in part supported by observations of quasar-driven 
outflows that find evidence for momentum and energy content 
comparable to that required by models (see, e.g.. [Tremonti et al.l 
20071; iReuland et al.ll2007l; I Arav et al.ll2007l;lNesvadba et alj|200a 
Prochaska & HennawfeOOaT but remains highly uncertain. 

In contrast, a separate class of models have been proposed to 
explain the BH-host correlations (or a subset) without feedback 
or self-regulation. An example is this: because the observed cor- 
relation between BH mass and host bulge stellar mass is approxi- 
mately linear (Mbh oc Mb u i gc ), if a fixed average fraction (~ 10~ 3 ) 
of the galactic baryon supply is able to reach the center and fuel the 
BH, this particular observed correlation can be explained. The ef- 
ficiency of angular momentum transport as well as the competition 
between gas accretion and st ar formation, explains a small mean 
constant of p roportionality dBurkert & Silkl 1200 ll ; lEscalal 120071 : 
iLi et al.ll2007h . Alternatively, stellar capture rates by the central 
BH based on the average central stellar density profile also lead 
to an expectation for the averag e BH accretion rate (see, e.g., 
iMiralda-Escude & Kollmeierll2005h . 

Although it is possible for these models to explain the 
mean correlation between BH mass and host mass, a more strin- 
gent constraint comes from the surprisingly small observed scat- 
ter in that correlation, which observations suggest is at most ~ 
0.25 — 0.3 dex (in terms of the constrained intrinsic scatter in 
Mbh versus Mb uig e , cr, or host potential depth/binding energy; see 
Tremaine et al. 2002; Marconi & Hunt 2003; Haring & Rix 2004; 



Graham & Driver 2007; Hopkins et al. 2007b; Aller & Richstone 



20071) . It has already been argued that this scatter (compared to 



that in other correlations) suggests that it must be galaxy proper- 
ties, rather than e.g. the properties of the larger d ark matter halo, 
that a re important for whatever sets the BH mass dWvifhe & Loebl 
2005). In addition, the small scatter puts strong constraints on 
mode ls where the mean re lation evolves very strongly with red- 
shift ((Robertson et al. 2006), although more moderate evolution as 
has been suggested from some observations is more easily recon- 
ciled with the low-redshift observations (see e.g. Peng et al] 20061 : 



iHopkins et al.ll2006d : ISal"viander et all2007l ; |Woo et alj|2006h . 

But the observed scatter also places strong constraints on 
whether or not some non-feedback mechanisms could be respon- 
sible for the observed correlations. It is not hard to imagine, after 
all, that — on average — ~ 10~ 3 of the galactic gas supply in most 
galaxies loses sufficient angular momentum to reach the BH radius 



of influence and ultimately accrete. What is hard to imagine is that 
it would always be this same fraction to within better than a factor 
~ 2, regardless of the huge diversity in observed galaxy gas frac- 
tions, gas and stellar mass density profiles, and kinematic states. 
On large scales where viscous forces are sub-dominant and gas dy- 
namics are largely gravitational, simulations of a galaxy merger or 
interaction suggest that even small changes to orbital parameters 
(say, changing relative disk inclinations, which are expected to vary 
widely, by a couple tens of degrees) lead to very large (nearly order- 
of-magnitude) changes in the amount of gas that loses such a large 
fraction of the angular momentum; from the galactic perspective, 
- 10~ 3 and ~ 10" " are both small, gravitationally negligible gas 
masses, so why should one be "picked out" so narrowly? 

At the opposite extreme, some models have proposed that 
feedback is so strong that the BH directly determines the stellar 
mass of the galaxy on a short timescale, as opposed to the feed- 
back models already mentioned in which the BH feedback pri- 
marily regulates its own growth. Such a scenario implicitly re- 
quires that the BH form in the background of a gas-dominated 
galaxy. The BH feedback is then assumed to either blow out all 
but the "desired" mass, which will become the final stellar mass 
appropriate for the observed correlations, or to directly induce 
star formation , causing th e formation of the appropriate mass in 



star tormation, causing the lormation ot tne appropriate mass in 
stars (see e.g. iKind l2003l 120051 ; iGranato et alj|2004l ; Isilkll2005l: 
iBeeelman & NatrJl2005T In these models, then, the BH 



mass is 

not a function of galaxy mass; rather, galaxy mass is a function of 
BH mass. As such, the scatter in such an inverse correlation is a 
strong constraint on these models. Moreover, because the galaxy 
mass formed at some distance from the BH depends explicitly on 
the BH mass, the simple expectation is that galaxies of a given BH 
mass should have similar stellar mass profiles near the BH where 
the galaxy can most directly "feel" the BH mass, but with increas- 
ing scatter at larger radii, where various effects will introduce scat- 
ter in "how much" of the BH feedback successfully affects gas at 
that initial radius. 

In this paper, we use observed galaxy stellar mass profiles and 
observations of gas masses, together with simulations of gas in- 
flows in galaxy interactions and mergers, to constrain the origin 
of the MeH-host galaxy correlations. We show that both observed 
galaxies and simulations exhibit a large scatter in the amount of gas 
that loses sufficient angular momentum to reach very small radii. 
Despite this order-of-magnitude scatter, the BH preserves (or sets) 
its mass such that it is correlated with the host on much larger scales 
to better than a factor of ~ 2. The fact that the total stellar mass on 
small scales has much larger scatter than on large scales is hard 
to reconcile with most existing models in which the BH is simply 
sensitive to the local mass supply. Indeed, this suggests that some 
process such as feedback should be present in order for the BH to 
be sensitive to global quantities such as the total mass of the galaxy 
(or total binding energy/central gas potential). The systematics in 
the radial distribution of scatter in the stellar mass also present a 
strong constraint on the most extreme feedback models, in which 
galaxy mass is set by BH mass, because these models must explain 
why the stellar mass at small radii is a less accurate tracer of the 
BH mass than at large radius. 



Even if the trigger for rapid BH growth is a galaxy merger, the back- 
ground potential at the galactic center relaxes sufficiently quickly that it 
is effectively fixed at the growth s tage of interest (see, e.g., lHopkins et all 
2007a, 2008di: lYounger et al.l2008l) . 



2 This is distinct from more common models of "radio" or "quiescent"- 
mode AGN feedback, in which AGN feedback regulates the final stellar 
mass of the galaxy over long ti mescales b y heating gas in a massive group 
or cluster halo (see e.g. lCroton et ah. 2006). 
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2 OBSERVATIONS: THE SCATTER IN THE ENCLOSED 
STELLAR MASS FRACTION IN ELLIPTICAL 
GALAXIES 

If the central BH mass were set by a fixed small fraction of the host 
mass reaching the central region, or if BH growth was determined 
by the central gas supply, we should expect that gas supply to be 
similarly well-correlated with the galaxy mass as is the BH mass. 
Although it is not possible to directly measure the total gas content 
that reached a given radius during BH growth, we can estimate this 
quantity and set a strict lower limit on it by measuring the stellar- 
mass inside a given radius in local spheroids. To lowest order, this 
will simply trace the total gas content that reached that radial scale. 
We discuss possible corrections to this below, but note that, strictly 
this sets a lower limit to the gas content that has entered a given 
galactic radius. 

Given a galaxy radial stellar surface density profile of the form 
E(i?), we integrate to determine the mass within R, M(< R). Since 
we are interested in comparing with BH mass, which follows a lin- 
ear relationship with total galaxy stellar mass (Mbh / M tot = /jbh ~ 
0.001 2 is assumed for all galaxies in the sample; see Harin g & Rixl 
2004), we divide out the total mass and obtain the quantity of in- 
terest, the mass fraction within R, f(<R)=M(< R) /M m ~ L(< 
R) /L to t (the last equality being for the observed surface brightness 
profile in any band where the mass-to-light ratio is a weak function 
of radius). Specifically 



/(<*)• 



M(<R) = J R X(R')2TTR'dR' 
M M ~ J °° £(#') 2-kR'AR' 



(1) 



The quantity /(< R) is at fixed physical size: systems with 
different effective radii will have significantly different /(< R) at 
the same R. Because the effective radius R e varies significantly as 
a function of total stellar mass, this would alone imply (even if 
galaxy profile shapes were identical) a very large variation in a 
population of galaxies in /(< R) at small R. Even at fixed stel- 
lar mass, the factor ~ 2 — 3 scatter in R e (M tot ) implies more than 
an order-of-magnitude scatter in f(< R) at small R if all galaxies 
have identical r'^-law fde V aucouleursl d 19481) profiles. Moreover, 
if the mass of the BH is a constant fraction of the mass of the host, 
then the BH radius of influence — the radius that should matter for 
the gas supply available for accretion — is a constant fraction of R e 
(since this radius is where the potential contributed by the BH and 
galaxy stellar mass are similar, i.e. GMbh/Rbh ~ GM tot /Re, hence 
Rbk ~ fiBHRe)- It is therefore more appropriate (and allows us to 
compare all objects on equal footing) to consider the mass fraction 
within a fraction of the effective radius, /(< R/R e ) (equivalently, 
to refer to all radii in units of the effective radius of the galaxy, 
rather than fixed physical units)(3 In short, for any fraction of the 
effective radius (or multiple of the BH radius of influence), we de- 
termine the fraction of the galactic gas supply that was available 
to turn into stars /(< R/R e ) by measuring the stellar mass fraction 
inside that R/R e . 

The quantity of interest is, however, not the median value of 
/(< R/R e ) (that is only a restatement of the typical profile shape 
of ellipticals), but how much scatter there is in /(< R/R e ) &IR/R e . 
Specifically, we consider a subsample of galaxies (usually at fixed 
stellar mass, to further marginalize over possible systematic dif- 
ferences, although we find this makes little difference), and de- 
termine /(< R/Re) by integrating the surface brightness measure- 
ments available, from the minimum resolution elements out to R e . 



By definition the mass fraction inside R e , /(< R/R e = 1) = 1 /2. 



For all the objects we consider, measurements extend to radii very 
much larger than R e . 

We then determine the scatter in log [/(< R/R e )], a[M(< 
R)/M tot ], in that sample at each R/R e , either by assuming the dis- 
tribution is lognormal and fitting it as such, or by taking the IPV 
value to reduce bias from outliers or skewness. We find it makes 
no difference. In terms of the simple dispersion (rather than a more 
complex fitted IPV value) this is 



i]T(log/,-<l°g/>) 2 , 



(2) 



with / = /(< R/Re) evaluated for all A' galaxies at the same R/R e , 
and (log/) the mean / at this R/R e . This procedure effectively 
gives a minimum scatter in the stellar mass fraction that has formed 
at R/Re or R/Rbh- Re-deriving the results that follow in terms of 
absolute physical radii (at fixed M 10t ) we find the same qualitative 
results with systematically more scatter at all radii, for the reasons 
given above. 

Figure Q] shows the results. We begin with a large sample 
of observed elliptical/ sphe roid surface bright ness profiles from 
iKorm endv et all ( 120091) and lLauer etail d2007l) . ~ 180 unique lo- 
cal ellipticals with nuclear HST observations and ground-based 
data at large radii (allowing accurate surface brightness profile 
measurements from ~ 10 pc to ~ 50kpc). Conversion to stellar 
mass profiles and comparison of profiles fro m different instruments 
and wa v ebands are disc ussed extensively in Hopk ins et all J2008b, 
I2009alldll2008l. r2009c); for our purposes the results are identical. 
We consider the scatter in mass fraction enclosed a[M(< R) /M tot ] 
as a function of R/R e in four bins of stellar mass, where there are 
sufficient numbers of observed objects (at least ~ 40 per bin) to 
obtain robust constraints. 

There is a clear mass-independent trend of scatter with R/R e , 
which we find can be reasonably well approximated as 



a[f{< R/Re)] ~ -0.28 log (R/Re 



(3) 



The scatter must go to zero, by definition, at R = R e because, by 
definition, all galaxies have exactly half their mass inside this ra- 
dius, but the rise towards smaller radii reflects a physical diversity 
of profile shapes. Repeating this exercise in terms of physical radii 
in small mass bins, a nearly identical trend is recovered, with a sys- 
tematically higher scatter at all radii by ~ 0.2 dex, and asymptoting 
to a constant ~ 0.2 dex scatter at radii corresponding to the typical 
R e at that mass and larger. 

Moreover, we have calibrated the expected results based on 
some simple experiments. For example, consider the case where 
we begin with a constant density profile over a number of fixed 
bins in R (analogous to the observed points) for some number of 
test cases, then independently perturb S(^) with a lognormal fixed 
scatter (say ~ 0.3 dex) at each R, and using the new density profiles 
for each system reconstruct Figure[T] We recover the "appropriate" 
answer - namely, that while the scatter does go to zero as R — > R e , 
the scatter at R <C Re converges to the input lognormal scatter as 
the number of mock systems is increased. 

Figure Q] shows that at R ~ 0. 1 R e , there is a factor 2 (0.3 dex) 
scatter in the interior stellar mass fraction, and that this grows to a 
factor ~ 4 (0.6dex) at R ~ 0.01 tf e and factor ~ 8 (0.8 dex) scatter 
near the BH radius of influence, R ~ 0.001 R e . In other words, by 
the time one is inside just 10% of R e , there is already more variation 
(galaxy-to-galaxy at fixed stellar mass) in the amount of mass that 
has "made it in" to smaller R than there is scatter in BH masses; this 
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Figure 1. Top Left: Scatter in the enclosed stellar mass fraction versus R/R e for ellipticals from the samples of Kormendy et al. ( 2009) and lLauer et alj (2007), 
in bins of stellar mass. Note a[M(< R)/M tot ] = at R = R e , by definition. At small radii, the scatter increases significantly. Top Right: Same, but for the most 
well-studied mass range (~ L»), and with the sample divided into nuclear "cusp" and "core" galaxies. Although the median central profile shapes are different, 
cr[M(< R)/M to t] is relatively independent of cusp/core status, suggesting no more than weak depend ence of a on subsequent ev olution and "scouring." Bottom 
Left: The "cusp" result compared with mass profiles of observed recent gas-rich merger remnants iRothberg & Joseph1 l2004) . The profiles in scatter agree, 
suggesting that the diversity in stellar mass at small scales is put in place by the recent star formation in gas inflows. Bottom Right: Comparison between the 
"cusp" ellipticals and the relaxed remnants of a diverse set of simulated gas-rich mergers (see ij3). In all cases, cr[M(< R)/M tot ] increases at small radii. 



difference only grows more and more pronounced as one moves to 
smaller radii. 

It is perhaps possible that some of this diversity on small scales 
owes to evolution subsequent to the formation and growth of the 
BH itself. For example, the difference between ellipticals with cen- 
tral "cusps" and those with nuclear "cores" (steep versus flat central 
profile slopes, respectively) is commonly attributed to "scouring," 
or the scattering of stars from an initially cuspy profile in three- 
body interactio ns by a binary BH af ter a "dry" or dissipationless 
galaxy merger (Bege lmanetaiT l980). To check this, we repeat our 
comparison in Figure [TJ (upper right panel) separately for galaxies 
classified as "cusp" and "core" systems (classifications for each ob- 
ject are given along with surface brightness profiles). To obtain the 
most robust statistics, we consider a bin in total stellar mass near 
~ L, (Mtot = 10 10 - 10" M Q ), where our sample includes ~ 30 
each of the cusp and core populations. At higher and lower masses, 
the comparison is similar, but the number of cusp and core ellipti- 
cals, respectively, drops rapidly. We find that there is no dramatic 
difference between the scatter within each of the two populations: 
the scatter in the amount of enclosed stellar mass increases dramat- 
ically at smaller radii for both cusp and core ellipticals. 

Alternatively, we can compare the results for galaxies that 



are known to be recent gas-rich merger remnants. Specifically, 
we repeat our analysis for the sample of confirmed gas-rich 
merg er remnants with near-inf rared surface brightness profiles 
from IRothberg & Jose ph (2004), al so discussed in the sam e con- 
text as the elliptical profiles above in Hopkins et aDH008b). These 
objects are relatively young, so have not had much opportunity to 
be affected by subsequent evolution, although they are sufficiently 
evolved s uch that they are dyna mically relaxed at R < R e . As dis- 
cussed in lHopkins et aD {2008b) we exclude all objects with obvi- 
ously unrelaxed features such as e.g. shells or tidal tails inside these 
radii, and check that the dynamical times at the radii of interest are 
shorter than the mean stellar population/secondary burst population 
ages. Although the resolution for these objects does not reach the 
extremely small radii of the Virgo elliptical HST sample, the sam- 
ple still overlaps over more than two orders of magnitude in R/R e . 
The lower left panel of FigureQ]shows that the scatter derived from 
this sample agrees well with the trend seen in cusp ellipticals of the 
same mass. We have experimented both with the entire merger rem- 
nant sample taken together and that sample binned by mass, and we 
find that our conclusions and the quantitative and qualitative shape 
of the scatter are unchanged. 

We emphasize that our discussion of the scatter in the enclosed 
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Figure 2. As Figure [T] but showing the logarithmic scatter in the stellar 
surface mass density S (R) (equivalently, mass in some annulus at R) rel- 
ative to some median profile shape for ellipticals of similar stellar mass 
(2(i? | M tot )), as a function of radius R relative to the effective radius of that 
median profile (R e (M tot )) (the median R e of spheroids of the given mass). 
In a narrow stellar mass bin, this is equivalent to the scatter in T,(R) at each 
R. Dividing the sample as in Figure[T]yields similar results. The scatter does 
not vanish at (R e (M t oi)), but it small - this reflects the scatter in the size- 
mass relation. The rise in scatter to small radius is similar to Figurefl]- the 
trends reflect genuine increased scatter in central properties, not an artifact 
of the fitting/quantities plotted. 



stellar mass with respect to the Mbh — host correlations implicitly 
assumes that all galaxies in the sample lie on the observationally es- 
tablished Mbh— host correlations, even though for many systems a 
BH mass has yet to be established. However we have examined the 
~ 10 cases in the adopted s amples that do have directly determined 
BH masses (compiled from Tremain e et al. ll2002l;lMarconi & Huntl 
12003k lHaring & Rixll2004l : lAller & Richstondl2007h and find that 
they are consistent with the above conclusions. The observed BH 
mass and radius of influence (calculated from the observed BH 
mass and velocity dispersion profile) agree with the simple esti- 
mates based on a constant BH to host mass ratio above. The scatter 
in e.g. the ratio MBH/M ga i(< R) exceeds an order of magnitude for 
R < R e , as predicted in Figure 1. For R > R e (M(< R) -> Mtot), 
on the other hand, the scatter in MBH/M ga i(< R) approaches the 
canonical ~ 0.3 dex. Bear in mind, however, that the number of 
such systems is small - rigorously re-fitting a BH-host mass rela- 
tion as a function of mass in different radii and estimating the in- 
ternal scatter will require larger overlapping samples with directly 
determined BH masses. 

As noted above, the definition adopted implicitly gaurantees 
the scatter vanishes at R = R e . In the interests of looking for a loca- 
tion where the scatter is "minimized," it is useful to define a similar 
quantity without this feature (even if the interpretation is somewhat 
less intuitive). For example, the scatter in surface stellar mass den- 
sity relative to some median profile of galaxies at the same mass. 
We show this in Figure[2] 

At a given stellar mass, we define the "median" profile as a 
Ide Vaucouleursl ([l948) r l//4 law with the same total stellar mass, 
but an effective radius equal to the median effective radius (R e ) of 
galaxies at that mass (from a quadratic fit to the R e — M tot relation 
in the sample). Knowing this (R e (Mtot)), the median profile is then 



where 64 and Z?4 are the appropriate normalization constants. We 
can then consider the ratio of the actual stellar mass surface den- 
sity T,(R) to E(if|M to t). Figure [2] plots the logarithmic scatter 
in this ratio at fixed value of R/ {R e (M tot )) (i.e. fixed R relative 
to the median R e at a given mass): i.e. Equation [2] where / = 
E(i?) /fj(R I Mtot) in some narrow interval (here 0. 1 dex intervals) in 
R/ (R e (M tot )). Because we are interested in the scatter, the median 
profile shape is implicitly normalized out; the results in Figure [2] 
are nearly unchanged if we assume a different Sersic profile with 
e.g. n — 2 — 8, or construct a non-parametric mean profile. 

The results are similar to Figure Q] the scatter in surface den- 
sity at a given radius (clearly related to the scatter in enclosed mass 
inside some radius) scales in a similar manner, reaching an order 
of magnitude at sufficiently small radii. Unlike FigureQ] the scatter 
does not have to go to zero at R — R e ; rather, at R = (R e (M tol ) ) , the 
scatter reflects th at in the size-mass relation of ellipticals (see e.g. 
IShen et alj|2003l) . Nevertheless, the scatter is much smaller around 
R e (and may even reach a minimum at these radii - rising again at 
/? 3> R e , reflecting the diversity in outer profile shapes where e.g. 
the presence or absence of extended envelopes is important). Like- 
wise, the conclusions are the same as Figure[T]if we split the sample 
by stellar mass or cusp/core status. 



3 THEORY: THE SCATTER IN THE ENCLOSED 
STELLAR MASS FRACTION IN SIMULATIONS 

We can gain some insight into how this scatter arises, and what 
it means relative to the amount of gas supply that has been avail- 
able to the BH, by comparing with the results of numerical sim- 
ulations of galaxy interactions and mergers that drive gas in- 
flows, and that form realistic elliptical galaxies. Speci fically, we 
consid e r the sample of hydrod ynami c simulations i n ICox et al.l 
(2006), iRobertson etal] d2006[). and lYounger et al.l J2008h a na- 
lyzed in detail in lHopkins et al.U2008bl. 12009 >J(j|2008all2009c1 rbl). 
These amount to several hundred unique simulations, spanning 
a wide range in progenitor galaxy masses, gas fractions, orbital 
parameters, progenitor structural properties (sizes, concentrations, 
bulge-to-disk ratios), and redshift. Most of the simulations are 
major (mass ratios ~1:3— 1:1), binary encounters, but they also 
include a series of minor mergers (from mass ratios ~1:20 to 
~1:3), as well as spheroid-spheroid "re-mergers" or "dry merg- 
ers" (i.e. mergers of the elliptical remnants of previous merger sim- 
ulations), mixed-morphology (spiral-ellip tical) mergers (see also 
Burk ert et alj|2008l : Ijohansson euiil 2009), multiple mergers, and 
rapid series of hierarchical mergers. Our results presented below 
are robust to these choices in the simulations. 

The numerical calculations usually include accretion and feed- 
back from supermassive black holes, as well as feedback from su- 
pernovae and stellar winds. However, we have performed parameter 
studies in these feedback prescriptions, and find that the structural 
properties of interest here are relatively insensitive to these effects 



E{R\M tot )=B 4 - 



Mu 



'Re(M tot )Y 



exp < 



{-K(^) 1/4 } (4) 



4 We determine the stellar surface mass density profile for each object from 
its surface brightness profile by assuming a constant stellar mass-to-light 
ratio as a function of radius and normalizing to the total stellar mass. The 
stellar masses are determined from the multi-band photometry (integrated 
until the light is converged), using the color-dependent mass-to-light ratio 
calibrations in Bell et al. (200 3j), assuming a " diet" ISarpeted 119551) IMF. 
Details are discussed in Hopkins et al. (2009a); changing the IMF, bands 
used, or allowing for a radius-dependent stellar mass-to-light ratio (based 
on the observed color gradients) makes little difference. 
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dCox et alJl2009l ; lHopkins et al.ll2007d.l2008bh . These calculations 
broadly reproduce the mean and individual profile shapes of objects 
in the observed samples, as a function of stellar mass, age, mor- 
phology, and galaxy type, as well as the fundamental plane scaling 
relations of spheroids. Since this sample reproduces the median and 
scatter in profile shapes, we have some confidence that it represents 
a reasonable proxy for the amount of gas inflow in real ellipticals. 

Figure Q] (bottom right panel) compares the scatter in interior 
mass fraction <r[M(< R)/M tot ] for this sample of simulations with 
the result for observed cusp ellipticals (blue dotted line). We con- 
sider several sub-samples of simulations to see how this scatter de- 
pends on input parameters. First, we define an "unrestricted" sam- 
ple that includes our entire suite of simulations. Second, we include 
a "restricted" sample that picks out a subset of nearly identical sim- 
ulations. These simulations have galaxies with similar gas content 
at the time of the final merger, identical stellar masses, include only 
equal-mass mergers, and are chosen from a narrow range in orbital 
parameters. The scatter in this case comes, therefore, only from es- 
sentially random object-to-object variance in the exact dynamics 
of the merger, and consequently in the exact amount of gas that is 
stripped of angular momentum at various radii. That the scatter in 
the resulting stellar mass fraction is already almost as much as that 
in the "unrestricted" sample demonstrates that most of the scatter 
in the full sample is a consequence of the random processes in a 
merger, rather than a result of the initial differences in gas fraction 
or orbital parameters. Thus, it is difficult to marginalize the scatter 
by invoking, for example, a narrow distribution in these parameters 
at a given BH mass. 

Third, we consider simulations with v ery strong stellar winds; 
these simulations are discussed in detail in lCox et"al] {2009). They 
include a mass loading factor several times the star formation rate 
and high wind velocities — effectively representing the m aximum 
wind feedback strength allowed b y observati o ns (se e e.g. [Martin 
ll999l ; lErbet"ai]|200d ; lirbl 120081) . ICox et alj j2009h demonstrate 
that yet stronger winds yield results which are physically incon- 
sistent with observations: galaxies "blow apart" any gas concentra- 
tion so efficiently that the existence of a starburst becomes impossi- 
ble. Thus, these "strong winds" simulations represent a reasonable 
physical upper limit to stellar feedback. The resulting dependence 
of scatter on radius is similar in the resulting galaxies, but with 
systematically higher normalization. 

Finally, we consider a subset of simulations with two dif- 
ferent specific choices for the sub-resolution prescription for the 
ISM gas (specifically the ISM equation of state: effectively the 
pressure support attributed to feedback loops in the ISM such 
as star formation, supernovae, stellar winds, radiation pressure , 
cosmic rays, and other sources; see e.g.lMcKee & Ostriker][l977l ; 
Springel & Hernquist 2003; Thompson et a lj2005b . This is param- 
eterized conveniently with the parameter g C o S : <?cos = being an 
isothermal equation of state for the ISM, the situation most unsta- 
ble to e.g. clumpin g and gravitational fragmenta tion, and q cos = 1 
representing the full lSpringel & Hernqu ist ( 2003) equation of state, 
which (over the density regime of interest) acts effectively as P oc 
p 1 ' 3-20 , leading to an ISM more resistant to clumping and grav- 
itational instability, with stronger pressure forces resisting infall 
and less stochastic scattering or torquing of gas clumps. Our re- 
sults are similar in either case because the details of the gas physics 
or stellar feedback prescriptions are largely secondary with respect 
to gravitational torque s, which dominate the dynamics in a merger 
dHopkins et"afll2009bh . 

Given these comparisons, we can see how the scatter in stel- 
lar mass content — which the simulations accurately reproduce for 



R/R e > 0.001 — relates to the scatter in the gas content "avail- 
able" to fuel the BH. Figure [3] shows this comparison. Specif- 
ically, we compare the correlation between BH mass and host 
spheroid stellar mass from the simulations (see |Pi Matteo et all 
l2005l : lRobertsonet"aIll2006h , with a small lognormal scatter of 
~ 0.25 dex, to the correlation between a number of measures of 
the nuclear gas and the bulge mass. First (upper right-hand panel), 
we consider the maximum gas mass at any time during the merger 
that is available inside of a nominal small radius (here lOOpc; the 
results scale with radius in the manner seen in Figure[TJ at smaller 
radii, however, our simulations begin to suffer from resolution lim- 
itations). We obtain a nearly identical result if we consider the to- 
tal gas content that enters or is supplied to a given radius (and of 
course, the total gas content that ultimately stays in a given ra- 
dius and is turned into stars is shown in Figure QJ. The scatter is 
clearly much larger than that between BH mass and host bulge 
mass, specifically in this case 0.72 dex. Second, we consider the 
gas mass in this radius ("available" for accretion - i.e. having made 
it to less than 100 pc inside the center of the starburst/gas inflow) at 
the time of the peak accretion rate/intrinsic luminosity of the BH - 
i.e. just before feedback from the BH begins to regulate growth of 
the system. Here (lower left-hand panel), the scatter is even larger, 
almost a full order of magnitude (0.97 dex). Finally, in the lower 
right-hand panel, we consider the maximum gas mass inside 1% 
of the effective radius (or fixed multiple of the BH radius of influ- 
ence), analogous to Figure [TQ The scatter is similar to that in the 
fixed 100 pc radius, and to that expected from FigureQ] 0.69 dex. 

Regardless of how we define the central gas content available 
for BH accretion, it seems that the scatter in the nuclear gas supply 
is always much larger than the scatter in BH mass at fixed host 
mass, in these simulations. Identical results are obtained in terms of 
other host properties like the velocity dispersion or binding energy. 

It is interesting to compare these predictions to alternative, 
more extreme feedback models, in which the BH mass and result- 
ing feedback set the total stellar mass itself by inducing star forma- 
tion and/or blowing out all but some desired amount of gas, which 
will then (later) form stars. This contrasts strongly with the simu- 
lations above, where the BH forms in a largely fixed background 
potential of stars, and self-regulates is own growth (with the mass 
blown out by the BH being a small correction to the total bulge 
mass). In the models where the BH induces the bulge formation, 
the mass formed at some R is an explicit function of Mbh or BH 
luminosity: M(< R) is a function of Mbh, rather than the other 
way around. As such, it is difficult to explain why the central re- 
gions of the galaxy near the BH radius of influence show a large 
scatter with respect to Mbh, while the mass formed at large radii, 
where the BH gravity is negligible, appears to reflect Mbh accu- 
rately. Quantitatively, we can construct analogous - albeit highly 
simplified - predictions for these models to compare with the stan- 
dard self-regulated BH growth scenario simulated above. Specif- 
ically, we can repeat the calculations used to derive the expected 
BH-host mass relation in the models, allowing for some flu ctuation 
in gas profile shapes. For example, we follow lKingU2003h and as- 
sume that BH feedback drives a compressive shock through a large 
gaseous body, inducing the formation of a bulge on the observed 
relation. Specifically, we assume that the gas initially fo l lows a uni- 
versal profile p(r) (in this case isothermal or Hernquist] dl990h . but 
it makes little difference for any p cx r^ 13 with j3 < 3, appropriate 



5 Here, we consider the true three-dimensional effective radius of the sys- 
tem. 
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Figure 3. Results from simulated merger remnants. Top Left: BH mass versus host bulge mass. The intrinsic scatter is shown (0.25 dex). Top Right: Maximum 
gas mass (at any time during the merger) within lOOpc of the BH, as a function of host bulge mass. A similar result is obtained if we consider the total gas 
mass that at some point enters the central lOOpc, or the mass in stars in this radius, as shown in Figure [T] The intrinsic scatter is three times larger. Bottom 
Left: Same, but for the gas mass inside lOOpc at the time when the BH is growing at its peak rate. Bottom Right: Same, but for the maximum gas mass 
measured within a fixed fraction of the remnant effective radius, rather than at fixed physical radius. In these simulations, the BH mass is more sensitive to 
global quantities, such as the stellar mass, that set the local potential depth/escape velocity, rather than the local gas content. 



for the radii r < R e that we are considering), and that the feedback 
drives a col d, infinitely thin shell which forms stars according to 
the observed Kennicutt ( 1998) relation; the details of the shell evo- 
lution are derived in iKind (2003, 2005). Now allow for some noise 
spectrum in object-to-object density profiles p(r) = p(r) (l + 5(r)) 
(S(r) is the amplitude of fluctuations in density per logarithmic in- 
terval in r, which we can parameterize |f| as 5{r) oc r 1 ) - this can 
represent clumping, fragmentation of shells, slightly different gas 
thermal states, inflow rates, or other factors. Integrating gives a pre- 
diction for the scatter in the log of the enclosed mass, as a func- 
tion of radius: very approximately u[f(< R/R e )] oc 7 + 0.5 logR. 
This is straightforward to understand: as the blastwave propagates 
away from the BH, the "knowledge" of the BH mass (accuracy with 
which the feedback strength at R reflects the true BH mass) exe- 



6 Equivalently, if there is some characteristic size scale rg oc r v for fluc- 
tuations in p and the density is a local gaussian random field, then 7 = 
(1 — r;)/2. Equal magnitude of random density fluctuations/deviations per 
logarithmic interval in r corresponds to 7 = or rj = 1 (rg oc r). Constant 
physical size-scales for characteristic fluctuations give 7=1/2. If the fluc- 
tuations seen in Figure [T] were put in as an initial condition in the density 
profile, it would be 7 = —0.3 (r] = 1.6). 



cutes a random walk. It is, however, opposite the observed trend 
(Equation [3} for most reasonable input noise spectra. This clearly 
represents a strong constraint on these models. But it may be pos- 
sible to reconcile them with the observations - for example, by in- 
voking the observed trend as part of the initial conditions, impos- 
ing an initial spectrum of density perturbations that already resem- 
bles what is observed^ This may be possible if, as demonstrated in 
the numerical simulations considered, gas is torqued and can move 
gravitationally in clumps and streams through an already largely 
formed stellar/dark matter background. However recall that in these 
models the background must (by construction) be primarily a self- 
gravitating gas-dominated system. It is unclear whether or not there 
actually exist equilibrium gas density profiles that can support (in 
an object-to-object sense) variations with as steep a spectrum as 
7 ~ — 1 (for the observed scatter, this would imply that the same 
median equilibrium large-scale profile could support variations of 



7 In detail following the calculation above, the initial perturbation spec- 
tram would have to be tuned as a function of radius such that 7 = —0.8 
(r? = 2.6); the fractional amplitude of local density perturbations reflecting 
the location within the galaxy and characteristically about two orders of 
magnitude larger at ~ 100 pc than at the effective radius. 
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an order of magnitude in the density inside ~ 100 pc, for example), 
f not so, other radius-dependent sources of scatter would need to be 
invoked to reconcile these models and the observations. 



4 CONCLUSIONS 

We have compared the scatter in the stellar mass in the central re- 
gions of observed spheroids, as a proxy for the amount of gas which 
was available at a given radius for possible accretion by the central 
super-massive BH. The scatter in this proxy for gas supply, from 
galaxy to galaxy, at fixed stellar mass, increases rapidly at small 
radii. It is nearly an order of magnitude around the BH radius of 
influence. We show that this is true for both observed nuclear cusp 
and core ellipticals, and for observed recent gas-rich merger rem- 
nants (see Fig.Q] ^2}. 

Despite this order-of-magnitude variation on small scales, the 
observed correlations between BH mass and global large-scale pa- 
rameters of host spheroids (e.g., total stellar mass, velocity dis- 
persion, binding energy) exhibit only factor ~ 2 scatter. In other 
words, the BH seems to "know more" about the host on large scales 
than on small scales. This is a natural expectation in feedback- 
regulated models: BHs stop growing once they reach a sufficient 
mass/accretion rate to unbind material near them, regardless of the 
total fuel supply. However, this finding is difficult to understand in 
the context of classes of models in the literature that do not invoke 
feedback and self-regulation, but instead posit that a fixed, small 
fraction, of the total host mass is incorporated into the central BH. 
If this were the case — if this gas supply set the BH mass — then 
the scatter in the nuclear gas supply would likely be directly related 
to the scatter in BH mass, and it should decrease or stay constant 
towards smaller radii. In fact, this is not observed. 

As a check, and in order to provide context, we compare the 
observed scatter with hydrodynamic merger simulations that in- 
clude feedback from growing BHs (see lower right-hand panel of 
Fig.0 Fig. [3] and |3). We find that the observed trends are broadly 
reproduced by the simulations. At various points in a galaxy's his- 
tory, mergers and interactions efficiently remove angular momen- 
tum from the gas, such that some of the gas (at some initial radii) 
effectively free-falls until shocking near the galactic center. As a 
consequence, the scatter in the amount of gas that reaches a given 
small radius is set primarily by gravitational physics, rather than 
the details of the gas microphysics or feedback prescriptions. We 
demonstrate explicitly that the scatter is also a general result of even 
very similar initial conditions: at small radii, the chaotic nature of 
mergers and interactions is important, and even strong restrictions 
in the initial gas fractions, redshifts, and orbits of interacting sys- 
tems does not significantly decrease the scatter in gas supply on 
scales comparable to the BH radius of influence. We also find that 
the simulations validate our assumption that the enclosed stellar 
mass is a faithful proxy for the supplied gas mass. Of course, it re- 
mains to be seen if this will be true as well for more "quiescent" 
fueling mechanisms (e.g. galactic bars and minor mergers, or stel- 
lar mass loss), that may be important or even dominant in the AGN 
population. 

These findings relate to another interesting property of 
spheroids: the absolute value of the central/peak stellar surface den- 
sity of spheroids is observed to be relatively constant (factor ~ 3 — 5 
scatter) independent of e.g. total stellar mass or effective radius 
of the galaxy (see e .g. lLauer et ail 120071 : [Kormendv et al.l [20091 : 
Hopk ins etal.l l2009a). which vary widely. If BH mass were set 
purely by local processes, the simplest expectation from this fact 



would be that all galaxies should correspondingly have the same 
BH mass (and that, as a consequence, there should be very large 
scatter between BH mass and total galaxy mass or velocity disper- 
sion). Instead, it seems to be these integral properties, rather than 
local density, that correlate best with BH mass. 

The constraints we derive on the run of observed scatter in 
stellar mass also appear constrain models at the opposite extreme. 
In these, feedback from the accreting BH is so strong that it de- 
termines the stellar mass of the galaxy, either by directly inducing 
star formation — essentially triggering bulge formation — or by 
blowing out a mass in gas much larger than the final stellar mass 
of the galaxy. In these models, galaxy mass is a function of BH 
mass, rather than the other way around. The scatter in the BH-host 
mass correlations is therefore a direct constraint on these models: 
in particular, not the absolute magnitude of the scatter, but the trend 
in scatter with radius is of interest. The stellar bulge mass formed 
in some shell at a given radius around the BH depends explicitly 
in these models on the BH mass (or luminosity), since this is what 
is triggering star formation or removing the gas supply that would 
otherwise make many more stars. Intuitively, it is difficult then to 
explain what is observed in FigureQ] despite the fact that the mass 
in each annulus should be determined by Mbh (if this class of mod- 
els were true), the central regions of the galaxy show more scatter 
with respect to Mbh than the outer regions of the galaxy. In par- 
ticular, the bulge mass at radii close to the BH radius of influence 
appears to have relatively little knowledge of the BH mass, exhibit- 
ing an order-of-magnitude scatter at given Mbh- In contrast, the 
mass formed at large radii, near the effective radius of the galaxy, 
here the BH's gravity is completely negligible appears to trace the 
BH mass accurately. There may be particular initial conditions or 
other controlling parameters that explain this, but it is clear that the 
observed scatter can strongly constrain such models. 

The systematics of the observed scatter in the BH-host rela- 
tions therefore — perhaps moreso than the normalization or slope 
of the relations — suggests that BHs may indeed self-regulate 
at a critical mass determined by global host galaxy properties in 
the manner predicted by feedback models. At the very least, non- 
feedback models must be revisited. It is insufficient to predict the 
normalization and slopes of the BH-host correlations. They must 
also provide predictions that reproduce the small observed scatter at 
large radii and its significant increase at small radii. Likewise, mod- 
els of very extreme feedback, in which the BH mass is not set by 
galaxy properties but rather, galaxy properties are set by BH mass, 
must explain how the central regions of the galaxy closest to the 
BH appear relatively insensitive to the BH mass. There is also the 
constraint from the BH "fundamental plane" - that it appears that 
neither Mbh — a nor Mbh — Mbuige is most fundamental, but rather 
a combination that traces bulge binding energy (~ (M bulgc( j 2 ) - 7 ) 
(all models should predict, for example, that Mbh correlates with 
a at fi xed Mbuige and vice ve rsa, as demonstrated in the observa- 
tions in Hopkins et alj2007bl) . Thus far, the class of self-regulating 
feedback models appear to be most successful at simultaneously 
explaining these observations, but more work is necessary. 

Finally, we note that the simulations described in S}3] indi- 
cate that the total stellar mass is a better predictor of Mbh than 
e.g. the stellar mass inside some small radius R. This prediction, 
coupled with the results in Figure Q] indicate that future observa- 
tions of the galaxies that comprise the observed Mbh— host rela- 
tions should also consider the radius-dependence of those relations 
(i.e. Mbuige — Mtot(< R) or Mbh — &(R) correlations, as a function 
of R). Comparison of such correlations, and identification of scales 
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where the scatter in the relations is minimized, can significantly 
constrain models of self-regulated BH growth. 
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